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' The determination of the numerical value for the strong coupling constant from tau decays in 

C . perturbative QCD is briefly described. Main emphasis is made on the use of the renormalization 

scheme invariance of the theory for reducing uncertainties related to truncation of the series in 
■ the strong coupling constant. The analysis of the convergence of the series is presented and the 

possibility of the asymptotic growth at the NNNLO is discussed. 

T— I ■ 
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• In this talk I briefly describe one aspect of r physics - the use of data on r-lepton hadronic 

. decays for extracting a numerical value of the strong coupling constant that is mainly based on 
the recent analysis presented in ref. W{. The accuracy of experimental data is steadily improving 

^ ■ that provides a good opportunity for high precision tests of the theoretical description [^Oll^. 

. The underlying theory of strong interactions - QCD - describes the observables related to the 

I physics of r-lepton hadronic decays within perturbation theory in the strong coupling constant, 

Q_i' i.e. as a series expansion. The key question of theoretical description is a pattern of convergence 



of the related series. On general basis, it is believed that the series is asymptotic, however, in 



D , some concrete applications it can appear as a convergent series. Thus, establishing the accuracy 

' with which the series, truncated in a given manner, represents the "full" result of perturbation 

theory calculation is a central question of using pQCD in low-energy phenomenology. The r-lepton 
^ . hadronic decays is a good place to investigate this problem since the theoretical description of the 

^ I process is very clean. Indeed, the basic quantity is the two-point correlator of hadronic currents 

in the Euclidean domain which has been calculated with a very high degree of accuracy within 
perturbation theory [SJ [3 |H] (for a review see |2| ) . Nonperturbative corrections to the correlator are 
known to be small and under control within the operator product expansion ^Uj- For comparison 
of experimental data with theoretical calculations one can choose inclusive observables that makes 
pQCD directly applicable to the description of the system ^1 ^1 ^1 El . These features 
explain the interest attracted by the r-lepton hadronic decays for extracting numerical values of 
the parameters of the standard model at lower energies. These values can then be compared with 
the results of high energy experiments that provides a powerful consistency check of QCD (and the 
standard model as a whole) for energies from M^- = 1.777 GeV till Mz = 91.187 GeV (e.g. |T7|). 

I discuss in some detail the procedure of extracting numerical values of as from r-data in 
perturbation theory. Since the numerical value of the expansion parameter is not small at the 
Mr scale the contribution of higher order terms in the perturbation theory series can be significant. 
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There are arguments that the accuracy of finite-order perturbation theory is aheady close to its 
asymptotic hmit ^H]. Therefore a kind of resummation is necessary. A powerful technique is 
1/Np expansion (e.g. ^]) that leads to renormalon-type resummation [201 El 1121 (SHI 1211 I21| 
that can be formulated as an integral of the running coupling constant in momentum space. This 
technique is often used in the analysis based on Dyson-Schwinger equations |26[ l27| . The result of 
the resummation under the integration sign depends strongly on the interpretation of higher order 
terms, e.g |28| I29 ( l3()]. The resummation of contributions related to the running of the coupling 
constant in QCD was used as a basis for the tau lepton decays analysis in contour improved 
perturbation theory ^1 . 

The central point to be discussed in the talk is an explicit use of renormalization group in- 
variance of the theory in order to analyze the r-lepton decay rate in the most optimal way within 
perturbation theory. Renormalization group invariance is a fundamental property of perturbation 
theory in quantum field theory which is (formally) related to the freedom in defining the subtraction 
procedure |32j . Its explicit compliance may make a numerical analysis more reliable and stable. 
Note that it can be viewed in a broader context as an optimal choice of the starting point for the 
perturbative expansion, i.e. the way of splitting the whole Hamiltonian into two parts: the leading 
one that is considered to be large and "perturbation" that is considered to be small and can be 
taken within perturbation theory. The influence of the proper choice of the initial approximation 
on the accuracy of the analysis can be clearly seen in model examples |33| or in the framework 
of effective theories |341 135| . In QCD the renormalization group invariance allows one to formally 
perform the numerical analysis in any renormalization scheme because all schemes are connected by 
a renormalization group transformation 36, 37 . In the finite-order perturbation theory approach 
this equivalence is only approximate and broken by higher order terms which introduce numerical 
differences into the results obtained in different renormalization schemes for the same quantities. 
Generally there many ways of using perturbation theory calculations. One is to find direct relations 
between physical observables which are renormalization group invariant. Then perturbation theory 
calculations are just a purely intermediate step for finding formal relations between observables 
(see, e.g. jl8l I38j ) and no numerical analysis for renormalization scheme noninvariant quantities 
is performed. Indeed, let the perturbation theory expressions for two observables Oi^2 in a given 
scheme have the form 



and is scheme- independent. The difference r2 — ri takes the same value for calculations in any 
scheme. Another way of using the result of perturbation theory calculations is to extract numeri- 
cal values for renormalization scheme noninvariant quantities (as the coupling constant in a fixed 
scheme) and then use it for predictions of different observables. The truncation of the perturbation 
theory series leads to numerical violations of renormalization scheme invariance and plays an es- 
sential role. In our simple example this means that the relations in eq. are treated as quadratic 
functions of in some fixed scheme and the accuracy of extraction of the coupling constant value 
(and prediction of other observables) depends drastically on the scheme used, i.e. on the numerical 
values of the coefficients ri^2- This happens because the accuracy is estimated in a heuristic way by 
considering the apparent convergence or the magnitude of the last term. There is an analogy with 
gauge invariance (or any other symmetry) : one can either work with invariant quantities or fix the 
gauge and compute in terms of noninvariant quantities but the final result for physical observables 



Oi = as + rial + O(a^), ©2 = + r2al + 0{al) . 



(1) 



Then the perturbation theory relation between observables Oi^2 reads 




(2) 
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(in a gauge invariant sector) will be the same. In practical application it is realized, for instance, 
as calculation in covariant or external field gauges. 

The concrete example of this talk is the extraction of a numerical value for the coupling constant 
which is not an immediate physical quantity. By convention the reference value of the coupling 
constant that is used to compare between different experiments is fixed to be the MS-scheme 
one. However, this does not necessarily mean that for its extraction from a given experiment 
the numerical analysis should be performed in the MS-scheme. It can be numerically accurate 
to analyze the system in its internal scheme and after finding numerical values for the internal 
parameters translate them into the MS-scheme using renormalization scheme transformation. This 
program is based on explicit renormalization scheme covariance of the theory (as calculations in 
different gauges are only possible (i.e. meaningful) for the gauge invariant sector of gauge theories). 
Note that the analogy with gauge is not quite complete: gauge invariance is exact order by order in 
perturbation theory while the renormalization scheme invariance is only exact for a 'full' quantity 
while for truncated series it is only accurate precise with accuracy the order of the value of the first 
omitted term. Note also that gauge invariant quantities could not always be put into the same order 
of perturbation series. The approach used below, in practice, means that any scheme is allowed for 
the analysis of a particular quantity while the transition between the schemes is considered to be 
exact. Therefore numerical values obtained in the MS-scheme directly and through renormalization 
group transformations can differ. In fact, in the spirit of the perturbation theory the internal 
scheme results are most reliable physically and are more stable numerically than the results of the 
standard analysis in the MS-scheme. The numerical values for the reference MS-scheme parameters 
are obtained by a renormalization group transformation from the numerical values found in the 
internal schemes. Renormalization group transformation is a quite formal operation and can be 
controlled numerically. For instance, the change of scale is a renormalization group transformation. 

We consider a simplest example that exhibits relevant features -the normalized r-lepton decay 
rate into nonstrange hadrons Hs=o that is given by 

RrS=o = ^^J^r ^f-^r^ = Nc\Vud\^SEw{l + Sp + 5ew + Snp) (3) 

where A^c = 3 is the number of colors. The first term in eq. © is the parton model result, the 
second term 6p represents perturbative QCD effects. V^d is the flavor mixing matrix element [1]. 
The factor Sew is an electroweak correction term [SH] and 5ew = 0.001 is an additive electroweak 
correction [30]. The nonperturbative corrections are rather small and consistent with zero; 6np = 
—0.003 lb 0.003 (see e.g. |13|'). For numerical estimates the factorization approximation for four- 
quark condensates is essential: it has been studied within x space sum rules and found to be under 
control j41j . Note that recently the problem of duality violation for two-point correlators has been 
discussed |42[ I43j . However, no established quantitative estimates of that violation are available 
yet. This problem can affect the numerical value of the coupling extracted from the analysis 
because of the numerical change of the quantity 5p extracted from eq. (jHl). However, no established 
quantitative estimates of that violation are available yet. We concentrate on the perturbative 
part of the decay rate and numerical uncertainties related to the renormalization scheme freedom 
of perturbation theory. In this respect new possible corrections do not qualitatively affect our 
analysis. The corrections due to duality violation are of a new nature and they can be added 
independently to eq. ©. They would only change the input numerical value for the 5p within our 
approach. 

The value for the decay rate RtS=q has been measured by the ALEPH 2^ and OPAL |3] 
collaborations with results very close to each other. Using (for definiteness only) the ALEPH 
data -^^5=0 ~ 3.492 it 0.016 one obtains from eq. © a numerical value for the main experimental 
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parameter of the discussed observable 

<5pP = 0.203 ± 0.007 . (4) 

Eq. @ is the main experimental input that is used in the analysis below. 

Theoretically one starts with the determination of the differential decay rate of the r lepton 
into a hadronic system H{s) with a total squared energy s 



2 / 2s 



ds ^ V M?) I M2 



Pis) 



that is given by the hadronic spectral density p{s) defined through the correlator of weak currents. 
For the {ud) current jj^ (x) = ^7/^(1 — 75)^ one finds in massless approximation for the light quarks 

i /(Tj7(x)jr+(0))e^''^dx = {q,q. - q\.)U^^%% U^^\q') = J ^ (5) 

with p{s) ~ Im U^'^'^is + iO), s = q'^. 

Integrating the function Il^'^'^{z) over a contour in the complex q"^ plane beyond the physical cut 
s > one finds that for particular weight functions some integrals of the hadronic spectral density 
p{s) can be reliably computed theoretically [Sllini- Indeed, due to Cauchy theorem one gets 



(p Il{z)dz = / p{s)ds . 
Jc Jcut 



Using the approximation n'^'^'^(2)|2gc' ~ n*^(2)|2gc' which is well justified sufficiently far from the 
physical cut one obtains 

/ U^''^{z)dz = f p{s)ds = i U'^\z)dz 

Jc Jcut Jc 

i.e. the integral over the hadronic spectrum can theoretically be evaluated. In practice, n*'^(z) is 
computed within OPE, i.e. the function Il^^{Q'^) with = _g2 -g calculable in pQCD far from the 
physical cut as a series in the running coupling constant asiQ"^) with power corrections. Further 
improvements on theory side can be made - for instance, instanton-induced contributions can be 
added. The lattice approximation for the evaluation of the correlator 11*^ (Q^) beyond perturbation 
theory can also be used [1^ . 

The total decay rate of the r lepton written in the form of an integral along the cut 



r(r ^ Hs=oi^) [ M Vi ^ 2s N 



2 

is precisely the quantity that one can reliably compute in pQCD. 

The basic object of the theoretical calculation is Adler's D-function which is given by the 
representation 

In the MS-scheme the perturbative expansion for the D-function is given by 

vr \ IT J \ IT J \ IT J 
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with (see e.g. [0]) 

299 ^, , , 58057 779 ^, , 75^, , 
^i = ^-9C(3), ^^ = ^-— C(3) + yC(5). (8) 

The notation as{Q) = as{Q)/7T for the standard MS-couphng constant normahzed at the scale 
fi = Q is used. Numerically one finds 

D{Q^) = 1 + a,(Q) + 1.6398a,(Q)2 + 6.3710a,(Q)3 + ksasiQf + Oi4{Q)) . (9) 

The coefficient is known only partly j47j . The particular numerical value of ~ 25 is obtained on 
the basis of geometric series approximation for the series © and is often used in the literature jl6l 
1481 149j . Below I keep this coefficient for illustrative purposes to see the potential influence of this 
term on the numerical value of the coupling constant extracted from r-data. 

In the MS-scheme the perturbative correction 6p is given by the perturbation theory expansion 

5f = as + 5.2023a^ + 26.366a^ + (78.003 + A;3)af + 0(af ) (10) 

where the MS-scheme coupling constant Og = vra^ is taken at the scale of the r-lepton mass 
fj, = Mr = 1.777 GeV. Usually one extracts a numerical value for agiMr) by treating the first 
three terms of the expression in eq. (|10)) as an exact function - the cubic polynomial, i.e. one solves 
the equation 

as + 5.2023a^ + 26.366a^ = 5p ^ . (11) 

The solution reads 

TTal\Mr) = al\Mr) = 0.3404 ± 0.0073e^p . (12) 

This is a standard method. The quoted error is due to the error in the input value of S^^. It is 
rather difficult to estimate the theoretical uncertainty of the procedure itself. The main problem 
is to estimate the quality of the approximation for the (asymptotic) series in eq. ()10() given by the 
cubic polynomial in eq. (fTT|) . 

As a criterion of the quality of the approximation one can use the pattern of convergence of the 
series (|10|) which is 

<5pP = 0.203 = 0.108 + 0.061 + 0.034 + . . . (13) 

One sees that the corrections provide a 100% change of the leading term. Another criterion is the 
order-by-order behavior of the extracted numerical value for the coupling constant. In consecutive 
orders of perturbation theory (LO - leading order, NLO - next-to-leading order, NNLO - next- 
next-to-leading order) one has 

afiMr)LO = 0.6377, af{Mr)NLO = 0.3882, af{Mr)NNLO = 0.3404. (14) 

One obtains a series for the numerical value of the coupling constant 

af{Mr)NNLO = 0.6377 - 0.2495 - 0.0478 -.. . (15) 

Up to the next-to-next-to-leading order result (NNLO) we can take a half of the last term as an 
estimate of the theoretical uncertainty. No rigorous justification can be given for such an assumption 
about the accuracy of the approximation without knowledge of the structure of the whole series: it 
is taken for definiteness. The theoretical uncertainty obtained in such a way - Aa'f{MT-)th = 0.0239 
- is much larger than the experimental uncertainty given in eq. p2() . This is a challenge for the 
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theory: the accuracy of theoretical formulas cannot compete with experimental precision at present. 
Assuming this theoretical uncertainty we have 

afiMr)NNLO = 0.3404 ± 0.0239t/, ± 0.0073exp. (16) 

Theory dominates the error even if the estimate for its precision ib0.0239t/i is not reliable (heuristic 
and only indicative). Thus the straightforward analysis in the MS-scheme is not stable numerically 
and the naive estimate of the theoretical uncertainty is large. 

To highlight the essence of the problem let me choose a different coupling constant as an 
expansion parameter that is obtained by the simple RG transformation - the change of scale 
of the coupling along the RG trajectory 1 GeV. It is often called scale transformation 

which is still a subgroup of the full renormalization group. In terms of 0^(1 GeV) one finds a series 
as(l)+2.615as(l)2 + 1.54as(l)^ = (5pP. The solution for the coupling constant is as(l GeV) = 0.453. 
The convergence pattern for the correction ^p'' is (5^'' = 0.203 = 0.144 + 0.054 + 0.005 and for 
the numerical value of the coupling constant 0^(1 GeV) = 0.453 = 0.638 — 0.177 — 0.008. Should 
one conclude that now the accuracy is much better? What would be an invariant criterion for the 
precision of theoretical predictions obtained from perturbation theory, i.e. finite number of terms of 
asymptotic series? Thus, one sees that the renormalization scheme dependence can strongly obscure 
the heuristic evaluation of the accuracy of theoretical formulae in the absence of any information 
on the structure of the whole perturbation theory series. 

The use of the MS-scheme is not obligatory for physical applications. The MS-scheme has a 
history of success for massless calculations where its results look natural and the corrections are 
usually small. This is not the strict rule, however, and there are cases (like gluonic correlators |5()| ) 
where corrections dramatically depend on the quantum numbers of the operators. In fact, the MS- 
scheme is rather artificial. It is simply defined by convention (let us be remindful of the evolution 
from the MS-scheme to the MS-scheme which had its origin only in technical convenience 
From technical point of view, in practical calculations of massless diagrams of the propagator type, 
another scheme - the G-scheme ~ is the most natural one It normalizes the basic quantity 
of the whole calculation within integration-by-parts technique ~ one loop massless scalar diagram 
- to unity j53j . /3-functions coincide in both schemes. It could have well happened that the G- 
scheme would be historically adopted as the reference scheme because corrections in this scheme are 
typically smaller than that in the MS-scheme. However, for the tau system the direct (standard) 
analysis in the G-scheme fails. 

Note that strictly speaking any scheme is suitable for a given perturbative calculation. However, 
it can lead to unusual (or even unacceptable) results in a numerical analysis. The only criterion 
for the choice of scheme at present is the heuristic requirement of fast explicit convergence: the 
terms of the series should decrease. Clearly this is a rather unreliable criterion. It does not provide 
strict quantitative constraints necessary for the level of precision usually claimed for the r-system 
analysis. 

One can however work within a different paradigm that is independent of the scheme in what the 
actual calculation has been performed. One just extracts a scale that any observable generates due 
to dimensional transmutation in perturbation theory and which is its internal scale. It is natural for 
a numerical analysis to determine this scale first. For comparison with other channels one can then 
transform the result into a MS-scheme or any other reference scheme using the renormalization 
group invariance. This last step is done only for comparison with other experiments (or just for 
convenience; the system itself can be well described in its internal scheme without any reference to 
the MS-scheme). 

The running of the coupling is one of the central features of QCD and very important for numer- 
ical analysis. A parameterization of the trajectory can be done at infinite Q by using dimensional 
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scale of QCD. The renormalization group equation 



,,^A^a{fi')=f3{a{fi')), a = - (17) 



is solved by the integral 



m (g) ^ + /*'(J^-^).« (18) 

where the indefinite integral $(a) is normalized as follows 

Here ^2(0-) and /3(a) denote the second order and full /? function, or as many terms as are available, 
given by 

/J2(a) = -a2(/?o + a/3i), Pia) = -a\po + P^a + P2a^ + Ps^^) + Oia^) , (20) 

a is a generic coupling constant. The four- loop /3- function coefficient /?3 is now known in the 
MS-scheme |S11 

140599 445 , , 

83 = + C 3) = 47.228 .. . (21) 

^ 4608 32 ^ ^ ^ ^ ^ 

The adjustment of the integration constant in eq. (|18|1 defines the standard QCD scale A^: the 
asymptotic expansion of the coupling constant at large momenta — > 00 reads 



The solution p8() of the renormalization group equation (|17p describes the evolution trajectory 
of the coupling constant. The evolution trajectory of the coupling constant given by the solution 
(|18|) of the renormalization group equation H17() is parametrized by the scale parameter A and the 
coefficients of the f3 function (3i with i > 2 (see e.g. J36 ). The evolution is invariant under the 
renormalization group transformation 



with the simultaneous change 



Po 1 left invariant and 



a{l + Kia + K2a'^ + K^a^ + . . .) (23) 
A^ ^ A2e-'^i//5o ^ (24) 



P2 P2- kIPo + K2P0 - KlPl 

Ps ^ P3 + 4kIPo + 2k3Po + kIPi-2ki{3k2Po + P2)- 

If this transformation was considered to be exact and the exact /3-function corresponding to the 
new charge was used then it would be just a change of variable in a differential equation (|17j) 
or the exact reparametrization of the trajectory and hence would lead to identical results. 
However, the renormalization group invariance of eq. (|18|) is violated in higher orders of the coupling 
constant because we consistently omit higher orders in the perturbation theory expressions for the 
/3-functions. This is the point where the finite-order perturbation theory approximation for the 
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respective /3-functions is made. This is the source for different numerical outputs of analyses in 
different schemes. 

Then the paradigm is to use the optimal charge and extract the coupling using exact RG 
technique. One introduces an effective charge a^- through the relation dp- = = arj'K jSHllSSlEni 
EH EHl and extracts the parameter A,- which is associated with a,- through eq. (|T8|l . This is just 
the internal scale associated with the physical observable Rr- The effective /3-function is given by 
the expression 

= -aliPrO + Priar + Pr2al + PrS^ + ■ ■ ■) (25) 

with f3ro = Po, Pri = Pi, and 

Pr2 = -12.3204 , pr3 = -182.719 + ^3 . (26) 

The extraction of the numerical value for the internal scale A,- is done from equation H18|) with 
a,-(M^) = 6'^^ . The coefficient P^-s does not enter the analysis. The parameter A^ = Aj^ is found 
according to eq. PH) . The MS coupling at ;U = Mr is obtained by solving eq. (|TH|) for ag^M-r) with 
regard to ln(M^/Ag) which is known if Ag is obtained; the /3-function is taken in the MS-scheme. 
For consistency reasons we only use the MS-scheme /3-function to three-loop order since the effective 
/^-function Pj- is only known up to the second order, cf. eq. (|26)) . A N^LO analysis is possible only 
if a definite value is chosen for ^3. Some estimates are given below. 

The procedure is based on renormalization group invariance and one can start from the ex- 
pression for the decay rate obtained in any scheme. The only perturbative objects present are the 
/3-functions: it is the converge of perturbation theory for /3-functions that determines the accuracy. 
It also highlights the limit of precision within this procedure: the expansion for Pt is believed to 
be asymptotic as any expansion in perturbation theory. The asymptotic expansion provides only 
limited accuracy for any given numerical value of the expansion parameter which cannot be further 
improved by including higher order terms. The expansion used is presumably rather close to its 
asymptotic limit 

pr{ar) = -al + Aar - 12.32040^ + al (^-182.719 + hz^^ + 0{al) (27) 

with ttr ~ 0.2 at the scale Mr- The convergence of the series depends crucially on the numerical 
value of ^3. If /c3 had a value where the asymptotic growth starts at third order then further 
improvement of the accuracy within finite-order perturbation theory is impossible. 

At every order of the analysis we use the whole information of the perturbation theory cal- 
culation. Especially, the appropriate coefficient of the /?T--function is present. In the standard 
method the coefficient P2 enters only at order 0(af) of the r-lepton decay rate expansion. The 
described procedure can be called the renormalization scheme invariant extraction method (RSI) 
that means that it is based on renormalization scheme invariance. It does not mean that renor- 
malization scheme uncertainty is completely eliminated but that it put under control in terms of 
explicit convergence of the effective P function. It is also clear since itself is not a physical 
object and is not renormalization scheme invariant. In this respect we extract the noninvariant 
parameter using invariance of the physics in order to perform the numerical analysis in the most 
suitable scheme. The output of the analysis is transformed into a numerical value for ag according 
to the renormalization group transformation rules that are treated exactly for the given order of 
the effective P function. For the coupling constant in the MS-scheme in NNLO it gives 

af^\Mr) = 0.3184 ± O.OOeOexp (28) 
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which is smaUer than the corresponding value obtained within the standard procedure eq. (|12j) . 
The key question is how to estimate the quality of this result? The parameter which is really 
extracted in consecutive orders of perturbation theory within the method is the scale A,-. Because 
of the relation 

A, = A^e-^-2°2^2/2/3o ^ o.3147A^ (29) 

one can look at A^ directly. One finds 

A^Ij^o = 595 MeV, A^Iat^o = 288 MeV, A^lArjvio = 349 MeV (30) 

or, representing the NNLO result as a formal series, 

A.|ArArLo = 595 -307 + 61 -... MeV. (31) 

Note that at leading order the scales (as well as charges) are equal in all schemes. Therefore 
the leading order result {As\lo = 595 MeV) is not representative, only indicative. Assuming the 
uncertainty of Ag to be given by the half of the last term of the series (|31|1 one has A^ = 349=b31 MeV 
which leads to the numerical value for the MS-scheme coupling constant 

as = 0.3184;0:0}^7 . (32) 

This result is obtained from eq. ()18() with three-loop /3-function. Taking the average we find 

a, = 0.3184 ±0.0159. (33) 

This is better than the theoretical error of the standard result eq. (|16|) . Still the theoretical error 
should be considered as a guess rather than a well-justified estimate of the uncertainty. It is 
affected by higher order terms, e.g. by the contribution. Clearly the estimate = 25 is rather 
speculative. Therefore it is more instructive to determine the range of k-^ which is safe for explicit 
convergence of perturbation theory. If the actual value of ^3 will be discovered in this range then 
perturbation theory is still valid and will give better accuracy in NNNLO. If not, the asymptotic 
growth of perturbation theory series is already reached and its accuracy cannot be improved. 

We require that the last term should be equal to the half of the previous one. In the standard 
way (eq. IW^ we have 

1(78 + k3)as\ < ^ 26.36 ^ 13 (34) 

which for = 0.1 gives 

- 208 < < 52 . (35) 

In the RSI way (eq. l27|) we have 



|(-182 + ^fc3)a,| < ^12.32 w6 (36) 

which for a,- = 0.2 gives 

33.8 < fcj < 47.1. (37) 

This range is much narrower than that in eq. H35|). The effective scheme method is much more 
sensitive to the structure of the series as can be seen from eq. (P7|) . The actual precision depends 
on the actual value chosen for ks and it is rather premature to speculate about numbers. 

Still we show the worst result (in the optimistic scenario that k^ lies in the safe range) that 
can be expected within the RSI approach. In the RSI approach with k^, = 47 we find the scale 
parameter in NNNLO 

A^Itv^tvlo = 334 MeV. (38) 



9 



With /c3 = 34 one has 



A. kiVTVLO = 367 MeV. 



(39) 



Taking the average we have 

As = 350 ± 17 MeV (40) 

which is the best possible estimate if we require that the perturbation theory series for the f3r- 
function still converges (according to our quantitative criterion of convergence) . That results in the 
numerical value for the MS-scheme coupling constant found with four-loop /3-function from eq. (|18|) 

0.3133 < as < 0.3314. (41) 

Therefore our conservative estimate of the theoretical error in the optimistic scenario for the con- 
vergence of perturbation theory series in NNNLO reads 

as = 0.322 ± 0.009 . (42) 

One should keep in mind that the theoretical uncertainty still depends on the criteria chosen. This 
particular estimate comes from the requirement of the explicit convergence of the effective beta 
function. 

The value of the coupling constant can be run to the scale Mz = 91.187 GeV with the RG For 
the standard method one finds 

af{Mz) = 0.1210 ± 0.0008e:,.p ± 0.0006^ ± O.OOOlb (43) 

where the subscript exp denotes the error originating from 6'^^. The errors with subscripts c,b 
arise from the uncertainty of the numerical values of the charm and bottom quark masses that 
enter the evolution analysis. The running to this reference scale is done with the four-loop /3- 
function in the MS-scheme |54| and three-loop matching conditions at the heavy quark (charm 
and bottom) thresholds (HHl- The threshold parameters related to heavy quark masses are = 
fhcific) = (1-35 ± 0.15) GeV and = fhb{nb) = (4.21 ± 0.11) GeV (e.g. pi) where ffiqiii) is the 
running mass of the heavy quark in the MS-scheme. 
The running of a^^^ {M-r) from eq. H33|) gives 

af^^{Mz) = 0.1184 ± 0.00074e^p ± 0.00053^ ± 0.00005b (44) 

The theoretical uncertainty comes mainly from the truncation of the perturbation theory series. 
Taking the result of the NNLO analysis eq. (|33]) one finds 

Aaf^{Mz)th = 0.0019 (45) 

In the most optimistic scenario with the NNNLO analysis eq. (|42jl one has 

af^\Mz)N-^LO = 0-119 ± 0.001 . (46) 

Note that for the Cabibbo suppressed decays that can also be analyzed within the perturbation 
theory |61M62ll^ the approach looks quite natural as there are additional functions to be studied: 
two functions related to perturbative corrections. Then three coefficient functions should be 
analyzed together that allows one to factor out the renormalization scheme dependence to large 
extent 

Now I briefly comment on the extraction of the strong coupling constant within resummed 
perturbation theory. There are basically two possibilities to resum perturbation theory series that 
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have recently been analyzed in some details jHS] • One is to extrapolate the running of the coupling 
into the infrared region using a natural cutoff provided by analytic continuation (so called vr^ 
terms j66l I67j). This approach relies on the pattern of continuation of the coupling constant to 
the infrared region. Another approach is based on the integration along the contour in Euclidean 
domain, it bypass the potentially nonperturbative region along safe areas in the complex momentum 
plane and, therefore, is quite perturbative in nature [14| . Numerically the results of these two 
approaches are different. In fact, they are also nonequivalent mathematically and the difference can 
be explicitly calculated: it happens to be nonexpandable in alg that means that it is not noticeable 
within perturbation theory. As for the numerical results one fits the theoretical expression for the 
decay rate in the contour improved approach to the experimental result 6^"^ eq. Q and find 

a'^\Mr) = 0.343 ± O.OOdexp (47) 

within the renormalization scheme invariant extraction method described above i.e. with the in- 
troduction of the effective charge first [^. This value differs from the finite-order perturbation 
theory result eq. which is expected. 

To conclude, the procedure of extracting the numerical value of the strong coupling constant 
from r-data is briefiy described with main emphasis on the particular features that are related to 
the renormalization scheme invariance of the theory. 

The present work is supported in part by the Russian Fund for Basic Research under contract 
02-01-00601 and by INTAS grant. 
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